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ABSTRACT 

We present Very Large Array observations at 7 millimeters wavelength that resolve the dust emission 
structure in the disk around the young star TW Hydrae at the scale of the ~4 AU ( (//IG) radius inner 
hole inferred from spectral energy distribution modeling. These high resolution data confirm directly 
the presence of an inner hole in the dust disk and reveal a high brightness ring that we associate 
with the directly illuminated inner edge of the disk. The clearing of the inner disk plausibly results 
from the dynamical effects of a giant planet in formation. In an appendix, we develop an analytical 
framework for the interpretation of visibility curves from power-law disk models with inner holes. 
Subject headings: circumstellar matter — planetary systems: protoplanetary disks — stars: individual 
(TW Hydrae) 



1. INTRODUCTION 

The TW Hya system is thought to be a close ana- 
log of the early Solar nebula. At a distance of 51±4 pc 
(jMamaiek 2005), it is the closest known classical T Tauri 
star, and a suite of observational studies have shown 
that TW Hya harbors a massive disk of gas and dust. 
Scattered light observations at optical and near-infrared 
wavelengths reveal a surface brightness profile consis- 
tent with a nearly face-on, optica lly thick, flared disk 
extending to ^200 AU in radius (Roberge et al.l 120051: 
Weinberger et all 120021 : iKrist et all 120001 : iTrilhng et all 



2001h . Observations at millimeter wavelengths have de 



tected thermal dust emission and a variety of molecu- 
lar species, i ncluding "CO, ^- ^ CO, CN, HON, HCO+ 
and DCO+ (IWeintraub et al.l 119891: IZuckerman et all 
19951: iKastner et all Il997t Ivan Dishoeck et all l2003t 



Wilner et al.H2003l: lOi et"all I2004D . The dust also dis- 

S'ays signatures of grain growth up to centimeter scales 
Vilner et al. 2005), and perhaps substantially larger 
sizes. 

Detailed models of the TW Hya spectral energy dis- 
tribution (SED) pr ovide constraints o n many aspects of 
the disk structure (jCalvet et al.ll2003 ). including the ra- 
dial dependence of outer disk surface density and tem- 
perature, and a clearing of the inner disk within ^4 AU 
radius. Resolved interferometric observations of millime- 
ter and submillimeter dust emission are in good agree- 
ment with the structure inferred from th e irradiated ac- 
cretion disks mode ls that match the SED (jOi et al.ll2004t 
I Wilner et a l."2GG0V though the resolution and sensitivity 
at these wavelengths have not been sufficient to address 
the presence of the inner hole. 

The inner hole is indicated by two features of the SED 
(jCalvet et al.ll2002f) : (1) a flux deflcit from -2-20^m, in- 
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dicative of low (dust) surface density in the inner disk, 
and (2) a flux excess at ~20-60 /im, thought to origi- 
nate from the truncated inner edge of the disk, directly 
illuminated by the star. Similar spectral features have 
been recognized in other T Tau ri star S EDs (e.g. CM 
Aurigae and DM Tauri; see Cal vet et al.l l2005) and may 
signify an important phase in the evolution of circum- 
stellar disks. One exciting possibility is that a discon- 
tinuity in the inner disk is a consequence of the per- 
turbative gravity field of a giant planet. Theories of 
planet-disk interaction predict the opening of gaps in a 
disk as a result of the formation of massive planets (e.g. 
Lin fc Papaloizo ul '19861 |Brvden et allll999( ). However, 
Boss fc York3 Tl993. 199^ show that the interpretation 
of infrared flux deficits as central clearings is not unique 
and reproduce SEDs of accreting disks around low-mass, 
pre-main sequence stars with a combination of opacity 
and geometry effects in the unresolved system. Spatially 
resolved observations of disk structure are required to 
confirm the inference from spectral deficits of inner disk 
clearing. 

To probe the disk morphology on size scale s commen- 
surate with the 4 AU transitional radius of .Calvet et al.l 
(|2002f ). we have used the Very Large Array^ to observe 
thermal dust emission from TW Hya at a wavelength of 
7 millimeters. These observations show clearly a deficit 
of dust emission in the inner disk consistent with the 
predicted hole. 

2. OBSERVATIONS 

We used the Very Large Array to observe TW Hya 
at 7 millimeters in the most extended (A) configuration. 
The observations used 23 VLA antennas (several were 
unavailable due to eVLA upgrades) that gave baseline 
lengths from 130 to 5200 kA. The observations were con- 
ducted for four hours per night on 10, 11 February 2006 
and 7 March 2006, from 7 and 11 UT (0 to 4 MST), 
during the late night when atmospheric phases on long 
baselines are most likely to be stable. Both circular po- 
larizations and two 50 MHz wide bands were used to 
obtain maximum continuum sensitivity. The calibrator 

® The National Radio Astronomy Observatory is a facility of 
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Fig. 1. — The TW Hya 7 millimeter continuum emission observed with the VLA (left ) compared to simulated images generated from the 
ICalvet et al.l {2003) model of an irradiated accretion disk, truncated at the 4 AU radius indicated by the SED and including the directly 
illuminated inner edge (center), or extending in to the 0.01 AU dust destruction radius (right). The contour levels are —2,2,3x0.23 mjy 
(the rms noise). In each panel, the ellipse in the lower left corner indicates the 0'.'29 X 0'.'09, PA 2.1° synthesized beam. The cross marks 
the derived position of the central star (see text). The center and right panels also show the disk models at full resolution, in a logarithmic 
grayscale to display the range of intensities in the fainter, outer regions of the disk and the bright, thin wall at the inner edge of the 
canonical model. 



J1037-295 was used to calibrate the complex gains, using 
an 80-second fast switching cycle with TW Hya. The cal- 
ibrator Jl 103-328, closer to TW Hya in the sky, was also 
included in a few minutes of fast switching each hour 
to test of the effectiveness of the phase transfer from 
J1037-295. The phase stability was good during the ob- 
servations of 11 February, worse on 10 February, and 
much worse on 7 March. Using the AIPS task SNFLG, 
we pruned the data with phase jumps of more than 70° 
between phase calibrator scans. This procedure passed 
about 80% of the data from the night of 11 February but 
substantially less from the other nights. Therefore, in 
the subsequent analysis we have used data only from 11 
February, the most stable night. The calibrator 3C286 
was used to set the absolute flux scale (adopting 1.45 Jy, 
from the AIPS routine SETJY), and we derived 1.95 Jy 
for J1G37-295. The uncertainty in the flux scale should 
be less than 10%. 

3. RESULTS 

3.1. 7 mm Image 

Figure [T] shows an image of the TW Hya 7 millimeter 
emission, where a Gaussian taper (FWHM 2000 kA) has 
been used in the imaging process to obtain an angular 
resolution matched to the surface brightness sensitivity. 
The rms noise level in this image is 0.23 mJy, and the 
peak flux of 0.88 mJy corresponds to a brightness tem- 
perature of 46 K. Because of the southern declination of 
the TW Hya and the VLA antenna geometry, the beam 
is elliptical and the resolution is higher east-west than 
north-south. Inspection of Figure [T] shows that the 7 mil- 
limeter emission is clearly not centrally peaked, as would 
be expected for a disk that extends continuously inwards 
towards the central star. Instead, the image exhibits a 
double-peaked morphology, consistent with a nearly face- 
on disk with a central hole observed with an elliptical 
beam. An image of the test calibrator J1103-328 made 
with the same parameters is point-like, as expected. 

3.2. Radially Averaged 7 mm Visibilities 

The central hole in the TW Hya emission may be iden- 
tified even more clearly in the visibility domain, free from 



the effects of the Fourier transform process and non- 
linear deconvolution. To better show the brightness dis- 
tribution, we averaged the visibility data in concentric 
an nuli of deproject ed (u,v) distance, TZuv, as described 
in lLav et all (|1997f l. We use t he TW Hya dis k position 
angle and inclination found bv lQi et al.l ()2004f ) from CO 
line imaging of the outer disk (-45° and 7°, respectively). 
These values may not be valid in the inner disk, if e.g. 
the disk warps in the interior. However, as long as the 
disk remains close to face-on at all radii, the deprojection 
correction is small and therefore insensitive to the exact 
values of these parameters. 

For each visibility, the coordinates were redefined in 
terms oi TZ — \Ju^ + , the distance from the origin 
of the (u,v) plane, and </> = arctan — PA), the polar 
angle from the major axis of the disk (defined by the 
position angle, measured east of north). Assum- 

ing circular symmetry and taking into account the disk 
inclination i, the deprojected (u,v) distances parallel to 
the major and minor axes of the disk, are da = 7?.sin0, 
df, = TZcoscf) cosi, respectively, and the deprojected 
(u,v) distance is TZuv ~ \/da + d^. 

An important parameter in the averaging process is the 
position of the star, or the center of the disk. To exam- 
ine the radial distribution of flux at the smallest scales 
permitted by the data, particularly in the east-west di- 
rection of highest resolution, we must know the phase 
center to within a fraction of the radius set by the reso- 
lution, i.e., the position of the star must be specified to 
within a few hundredths of an arcsecond, which is better 
than the absolute astrometric accuracy of the data (typ- 
ical positional accuracies in the A array are ~ O'.'l due 
to baseline uncertainties and uncorrected tropospheric 
phase fluctuations; see, e.g. the 2004 VLA Observational 
Status Summary). To the extent that the disk is sym- 
metric, the process of deprojecting and averaging at the 
correct star position will minimize the scatter within each 
deprojected radial bin and bring the average of the imag- 
inary parts of the visibilities (the average phase) to zero. 
Therefore, we chose the star position to be that which 
minimized the absolute value of the mean of the imag- 
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Fig. 2. — TW Hya real and imaginary 7 millimeter visibilities, 
in 430-kA bins of deprojected u-v distance from the center of the 
disk. Error bars represent the standard error of the mean in each 
bin. Note that bins represented by open and filled circles are not 
independent: the bins are overlapping, i.e., the domain of each 
open circle extends to the neighboring filled circles. The calculated 
visibility functions for the best fit irradiated accretion disk model is 
indicated by the heavy solid line and is composed of contributions 
from an outer disk (solid line) with an inner hole of radius 4.5 AU 
and a bright, thin wall (dotted line). For comparison, the visibility 
function for a disk model that extends in to the dust destruction 
radius at ~0.01 AU is also shown (dashed line). To the extent that 
the disk is radially symmetric, the imaginary part of the visibility 
(the average phase) should be zero for all models. 

inary visibility bins. This position is indicated by the 
cross in the left panel of Figure [T] 

Figure [2] shows the annularly averaged visibility am- 
plitude as a function of TZ^v The width of each bin is 
430 kA, chosen to be narrow enough to sample the shape 
of the visibility function and also wide enough to have 
sufficient signal-to-noise ratio. Although the visibility 
data are still noisy when divided up in this way, it is 
evident that the visibility function passes through a null 
near TZuv of ~1000 kA, indicative of a sharp edge in the 
emission. 

4. DISCUSSION 

The most striking feature of the new 7 millimeter ob- 
servations is the central depression in the image, which 
is also indicated by the presence of the null in the depro- 
jected visibility function. We identify this feature with a 
clearing of the inner dust disk. A continuous disk that ex- 
tends inward to the dust destruction radius at ~0.01 AU 
would show sharply centrally peaked emission and would 
not show a null at the observed baselines. Figure [T] com- 
pares the 7 millimeter image (left panel) with an im- 
age generated for a model with a continuous disk (right 
panel). Figure [5] also shows the visibility function of con- 
tinuous disk (dashed line); this is clearly not compatible 
with the 7 mm observations (or the infrared SED ). 

Independent of detailed modeling, the size scale of this 
inner hole can be estimated simply from the separation 



of the peaks in the image. These peaks are sensitive 
primarily to inner edge of the disk, where the 7 millime- 
ter brightness is highest. The separation of the peaks is 
0'.'14, or ~ 7 AU. This separation is slightly smaller 
than the diameter of the inner hole, since the bright emis- 
sion from the inner edge to the north and to the south of 
the star, combined with the lower angular resolution in 
the north-south direction, tend to draw the image peaks 
together. This effect is evident in the right panel of Fig- 
ure [U in which the contoured peaks can be seen to lie 
interior to the inner edge of the disk. The null in the vis- 
ibility function provides a corroborating estimate of the 
size of the inner hole in the disk. In appendix [A} we show 
how the angular scale of the null in the visibility function 
of a power law disk depends on the density and temper- 
ature power law indices and the radius of the inner hole. 
If we assume that the outer disk emission contributes lit- 
tle on these long baselines, and that the total emission 
is dominated by a bright, thin ring associated with the 
inner edge of the disk, as discussed in i j4. 1.11 below, then 
we can estimate the radius of the inner hole hole with 
equation lAllI A linear fit to the binned visibilities gives 
a null position of 930±60 kA and implies an inner hole ra- 
dius 4.3±0.3 AU. Thus the resolved dust emission shows 
an inner hole in the disk at a size scale very similar to 
that inferred from SED modeling. 

4.1. Comparison with Disk Models 
4.1.1. Power Law Disk Models 

The dust emission from the outer disk of TW Hya 
has been shown to be a good match to the structure 
of an irradiated accretion disk model, approximated by 
power laws in temperature and surface density with in- 
dices q ~ 0.5 and p ~ 1.0, respectively, over a wide 
range of radii. However, an extrapolation of this very 
simple power law model to an inner disk truncated at 
~ 4 AU radius is not compatible with the new long- 
baseline 7 millimeter data. For these power laws, the 
null in the visibility function of the 4 AU hole should ap- 
pear at ^ 500 kA, which is a factor of two smaller than 
observed (equation IA9p . In order for a 4 AU hole to be 
consistent with a ~1000 kA null, the power law model 
requires a much steeper emission gradient, with the sum 
of the radial power law indices p-\-q approaching 5. Such 
steep power laws are inconsistent with observations of 
the outer disk. In addition, this power law model fails to 
reproduce the flux in the image peaks by nearly an order 
of magnitude. Another possibility is that the radius of 
the hole is smaller than that predicted by the SED. A 
power law disk with p + q ~ 1.5 and radius 2 AU does re- 
produce the 1000 kA null and observed peak separation; 
however, this model still fails to reproduce the observed 
peak flux by more than a factor of five. These discrep- 
ancies indicate that a more complex model of the disk is 
needed. In general, reproducing the observed peak flux in 
the new high resolution observations requires a greater 
concentration of material at the inner edge of the disk 
than that of power-law disk models. 

A natural modification to the power law disk model 
with an inner hole is the addition of a bright, thin, inner 
edge, or "wall" component. This "wall" component cor- 
responds to the frontally illuminated i nner edge of the 
disk in the calculations of lCalvet et al.l ()2002[ ). who show 
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that a small range of temperatures is required to repro- 
duce the narrow spectral width of the mid-infrared ex- 
cess. The presence of this additional compact component 
to the model shifts the angular scale of the null in the vis- 
ibility function to larger TZuv and raises the flux at long 
baselines. In this composite model, for a given power law 
description, the angular scale of the null in the composite 
"disk-l-waU" visibility function depends on (1) the radius 
of the inner hole, and (2) the relative brightness of the 
disk and wall. The effect of the second dependency is 
to move the angular scale of the null between the limit- 
ing positions from the "disk" alone and from the " wall" 
alone. In fact, the position of the null in the 7 millime- 
ter data is close to that expected from an infinitesimally 
thin ring of AU radius fequation lAll[) . A bright, thin 
ring also reproduces the ~ 0'.'14 separation of the image 
peaks. Thus it appears that the high resolution obser- 
vations show primarily the directly illuminated wall at 
the inner edge of the disk. At these long baselines, the 
extended emission from the outer disk that dominates at 
larger size scales is weak and effectively not detected. 

4.1.2. Irradiated Accretion Disk Model 



The irradiated accretion disk model of ICalvet et al.l 
(|2002t l provides a more realistic description of disk struc- 
ture than a simple power-law model. To compare the 
data to this more sophisticated disk model, we simulate 
numerically the expected emission at 7 millimeters, in- 
cluding the detailed visibility sampling of the Very Large 
Array observations. 

For the frontally illuminated component of the in- 
ner edge of the disk, we follow the prescription of 
ICalvet et al.l (|2002f ). adopting a shape given by 



zoexp[{R~ Ro)/AR] 

and temperature 

1/2 



Tphot{R) ~ 



R 



(1) 



(2) 



where R is the radius from the star, Zg is the height 
of the wall, subscript refers to the boundary be- 
tween the wall and the outer disk, = cos 6*0 and 
^0 = Tr/2 — ia.n~^{dzs/dR), and we assume that at 7 mil- 
limeters the radial brightness tracks the temperature. To 
set the absolute flux of this component, we normalized 
the intensity distribution to match the peak flux of the 
image and the total flux of the di sk determined from 
previous, lower resolution imaging (iWilner et al.l[200l . 
Since Ai? is not well constrained by our data, except that 
the wall must be narrow compared to the size of the hole 
(Ai? < K) and the resolution of the data (Ai? < 0'.'09, 
or 4.6 AU at 5 1 pc) , we use the value Ai? = 0.5 AU which 
ICalvet et a l.l (J2002I) find to be consistent with the shape 
of the mid-IR SED. We choose Rq to be 4.5 AU so that 
the peak of the wall emission occurs near 4.0 AU, causing 
the null of the model to match the position indicated by 
the data. For the dust mass opacity, we use the power 
law form k = ^o(i^)^ where kq = 1.8, vq = 1.0 x 10^^, 
and p = 0.8 ()D'Alessio et al.ll200H ). 

We u sed the Monte Carlo radiative transfer code RA- 
TRAN (jHogerheiide fc van der Tak.2000 ) to calculate a 
sky-projected image from the model continuum emission. 



with frequency and bandwidth appropriate for the obser- 
vations, and the miriad task uvmodel to simulate the ob- 
servations, using the appropriate antenna positions and 
visibility weights. Figure [2] shows the visibility function 
from this model (heavy solid line), which matches well 
the observations. The light solid line shows the contri- 
bution from the outer disk, which accounts for only a 
small fraction of the flux observed at long baselines, and 
the dotted curve shows contribution from the thin wall 
component, which dominates the outer disk at TZuv be- 
yond ~ 100 kA. The long baseline data are sensitive to 
the total flux and thin width of the wall, and are com- 
patible with the assumptions of the wall structure used 
in the previous SED modeling. In the imag e domain, we 
have c ompared the 7 millimeter data to the'Calve t et al.l 
l|200l model using RATRAN to generate images of a 
disk with an inner hole of radius 4.5 AU and a bright 
wall of width 0.5 A U and total flux 1.7 mJy (Figure [H 
center panel), which [Cal vet et all ()2002[ ) find to be con- 
sistent with the mid-infrared SED. For comparison, we 
also generate an image of a continuous disk extending 
inward to the dust destruction radius at 0.01 AU (Fig- 
ure [l] right panel). The data are inconsistent with this 
continuous d i sk mo del, while the SED-based model of 
ICalvet et all (|2002D with an inner hole and bright wall 
matches the observed structure well. 

In short, the 7 millimeter data suggest an outer dust 
disk extending inward to ^ 4.5 AU, where there is wall 
of radial extent ~0.5 AU that is much brighter than the 
surrounding disk on account of direct illumination by the 
star. Interior to the wall, a sharp transition occurs to a 
region of lower dust surface density and correspondingly 
weak or absent 7 millimeter emission. 

4.2. Disk Clearing 

The resolved observations bolster the inference from 
SED models that the TW Hya disk has an inner 
hole of much reduced dust column density of radius 
^ 4 AU. Theories of disk-planet interaction have long 
predicted the opening of gaps in circumstellar disks as 
a consequence of the formation of giant planets (e.g. 
iLin fc Papaloizoull986tlBrvden et allll999[ ). and numeri- 
cal simulations of such interactions produce these gaps in 
a consistent way acros s difltering architectures a nd com- 
putational algorithms fde V al-Borro et ani2006[ ). Recent 
work by Varnicrc ct_al. (200^ has shown that in turbu- 
lent disks with a-viscosity, inner holes are in fact more 
likely than gaps as a consequence of the formation of 
Jupiter-mass planets around solar mass stars, as the pro- 
cess of outward angular momentum transfer mediated by 
spiral density waves can cause clearing of the inner disk 
on time scales an order of magnitude shorter than the 
viscous timescale. 

While the inner disk is largely cleared, it is not entirely 
devoid of gas and dust. Optical emission lines indicate 
gas accretion (Muzcrollc et al. 2000 ), alb e it at the low 
rate of 4 x 10'^° Mq yr'^. IRettig et all (|2004[) detect 
'--^ 6 X 10^^ g of warm CO at distances of 0.5-1 AU from 
the star, and iHerczeg et al] ()2004f ) infer the presence of 
warm II2 within 2 AU of the star by modeling the HST 
and FUSE H 2 emission spectrum. The 10 /i m silicate 
feature fSitko et all 120001 : lUchida et al.ll200l indicates 
that there must be at least a few tenths of a lunar mass of 
dust grains present in the otherwise largely cleared inner 
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disk (jCalvet et al.ll2002D . T he region within ~0 .3 AU has 
been spatially resolved by lEisner et alJ ()2006f ) at 2 
using the Keck interferometer, detecting emission from 
optically thin, submicron-sized dust populating the in- 
ner regions. This small amount of material may be the 
result of a restricted flow through the transition region 
at ~ 4 AU from the massive reservoir of the outer disk. 
A population of micron-sized gr ains in the inner disk is 
consis tent with the predictions of lAlexander fc Armitaed 
(|2006) , who show that planet-induced gaps tend to filter 
out and confine large grains at the gap edge while allow- 
ing small grains to migrate across the gap with accreting 
gas and populate the inner disk. Our data are consistent 
with an inner region devoid of 7mm dust emission; how- 
ever, the degree of clearing is difficult to constrain due to 
both the low signal-to-noise ratio at long baselines and 
the dependence on the wall model adopted. 

The dynamical effect of a planet is not the only pos- 
sible explanation for the observed central flux deficit at 
millimeter wavelengths. Photoevaporation of gas by ul- 
traviolet radi ation has also been invoked to explain inner 
disk clearing (jClarke et al.ll2001h . if the evaporation rate 
at the gravitatio nal radius domin a tes th e accretion rate. 
As discussed by [Alexander et al.l ()2006[ ). TW Hya is at 
best marginally consistent with a photoevaporation sce- 
nario, since the outer disk mass is larger than predicted 
for photoevaporation, and the observed accretion should 
occur only for the brief period while the inner disk is 
draining onto the star. In contrast with the case of a 
planet-induced gap, photoevaporation should also clear 
the inner disk of all but the largest grains, as noted by 
[Alexander fc Armitaed ()2006D . Other mechanisms, such 



as photophoretic effects (jKrauss fc Wurml[200l . could 
also aid in clearing the inner disk region, if gas densities 
are high enough. 

5. CONCLUSIONS 

We present new, spatially resolved observations of the 
TW Hya disk at 7 millimeters that provide direct evi- 
dence for a sharp transition in dust surface density at 
~ 4 AU radius, a feature cons istent with the inne r hole 
inferred from SED modeling bv lCalvet et"aLl (|2002D . The 
interpretation of the mid-infrared flux deficit as a central 
clearing of material is robust. The TW Hya system is ide- 
ally suited to future observations that will be able to dis- 
tinguish between the various scenarios invoked to explain 
central clearing in the disks around young stars. Signa- 
tures diagnostic of planets in formation, in particular spi- 
ral density waves in the disk and thermal emission from 
circumplanetary dust, should be within the detection ca- 
pabilities of the Atacama Large Millimeter Array operat- 
ing at the shortest submillimeter wavelengths (300 /im) 
and at the longest baseline s (> 10 km) to achieve t he nec- 
essary ang ular resolution (jWolf fc D'Angelol[2f)n5h . Such 
observations will benefit from the close proximity of TW 
Hya, the nearly face-on viewing geometry, and the size 
scale of the inner hole now confirmed by direct observa- 
tion at millimeter wavelengths. 
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APPENDIX 

A. PROTOPLANETARY DISK VISIBILITY FUNCTIONS 

We analyze a simplified disk model parameterized by power law distributions in surface density and temperature 
and including a central hole, to illustrate how the model parameters affect the shape of the visibility function. We also 
discuss the visibility function of a model consisting of a thin ring. For both cases, we provide analytical expressions 
for calculating the position of the null in the (deprojected) visibility function, which is an easily observed feature. 



A.l. Power-Law Disk with a Central Hole 

For a flat, optically thin disk described by power-law distributions in temperature and surface density, assumed to 
be radiating in the Rayleigh-Jeans limit and viewed face-on, the intensity of radiation as a function of radial angular 
scale 6 from the center of the disk is 



I AO) = 



Jo 



■^0 



-V 



The visibility as a function of (u,v) distance, 7?., is given by the Fourier transform of the intensity: 

where Jq is a zeroeth order Bessel function. The integral term in the above expression can be evaluated as 

/)2-(p+9)p _ p+£\ 



where 



1F2 



{l-^},{l,2-^},-.W 



(Al) 



(A2) 



(A3) 



(A4) 



and xFy{a.,h, z) is a generalized hypergeometric function of one variable. For a disk with an inner hole, the limits of 
integration da and 9b are the inner and outer angular radii of the disk, respectively. 
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Fig. 3. — The dependence of TZ on disk temperature and surface density power law indices (p + q), source distance (Dsource), and inner 
hole size (Rhole)- 

For characterizing the visibility function on angular scales between 9a and 6*^, we can make the approximation 
Tl^O^ <C 1 and 7i?0^ ^ 1. In these limits, the hypergeometric function attains manageable analytic forms. 
For 2; <g; 1, appropriate for the inner disk of radius Oa, the function x-Fj/(a, b, z) has the following series expansion: 



lim ^Fy{a,h,z) 



PH{c,k) = {c + k)\/d (A5) 



and in this limit the quantity J- goes to 



2- {p + q) ne^^o 



with the limit approximated by the sum in IA5I above. 

In the limit of 2; 2> 1, appropriate for the outer radius Ti?6^, the quantity T approaches the analytical form 



(A6) 



lim TiOb) = — ^ A7 



The total visibility function (jA21 with integral IA3p is then the difference between a smooth power-law disk without a 
hole {T{9b), equation I A7p and the contribution of the evacuated inner region {!F{da), equation IA6I with limit IA5| ). 



A. 1.1. Position of the Null 

For a disk with an inner hole, the position of a null in the visibility function is an easily observed quantity. Here we 
show how the angular scale of the first null depends on the disk model parameters. Substituting expansion IA5I and 
equation IA7I into IA3I and setting the result equal to zero, we obtain the following expression which can be solved for 
the position of the first null: 

r (1 - ^) (1 - „ .2-ip+q) 2 - (p + q) _2-r>2 n2 , \ (Ao\ 
Y^P+q^ = i^T^nuuOa) [I - T^nulA + ' ' ' j (A8) 

Since TZ and 9a always appear in tandem in this expression, TZnuii oc j-. 

The dependence of the null position on the power law indices p and q is illustrated in Figure [3l which shows the null 
position as a function of p + g for a fixed inner disk radius. Several orders of the power series expansion (jASP are shown. 
For typical values of the power law indices, the null position shifts monotonically to longer baselines as p + <z increases, 
exhibiting an essentially linear relationship in the vicinity of p + q—2. As p + q increases, the temperature and surface 
density distributions (and therefore the intensity) become more sharply centrally peaked, and so the position of the 
null, which is effectively the angular scale on which the inner disk contribution to the visibility equals that of the outer 
disk, moves to smaller and smaller angular scales (i.e., larger TV). 

These approximations begin to lose validity longwards of the vicinity of the first null, which occurs at TZOa < 0.3 for 
allp + q < 3. The series expansion quickly diverges past TZ9a = 1. However, it should be also noted that this constraint 
places no limit on the size of the hole that can be investigated by this method, i.e., for any disk with a central hole 
there will always be at least one null shortward of TZ9a—l, and so this method is robust for any case in which 9f, is large 
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Fig. 4. — The visibility function for the TW Hya disk, based on the power-law model with an inner hole in which expansion lASl is carried 
to th ird order (solid line). The visibility function for a thin ring interior to the inner disk edge is also shown (dotted line; see equation 
RTOt . 

compared to 9a- It is also valid for an inclined (symmetric) disk, as long as the deprojection is handled appropriately, 
as in W6.2\ 

For ease of use, it is possible to approximate the {jp + q) dependence by a linear fit of the curve in the region 
1 < p + g < 3, which results in the following formula for the position of the null, good to within 4%: 



Ti-nullikX) 



fl AU\ f D 

V Rhoh 



100 pc 



[2618+ 1059 X {p + q)] 



(A9) 



A.2. Thin Wall 



A thin wall (Ai? ^ R) can be described by a ring of constant brightness at a distance Rhoie from the star, with a 
visibility function which is a zeroeth order Bessel function: 



/•oo 

V{n) = 27r / I^{e)Ja{2Tr9n)0d0 
Jq 

= 2Tr9iioleIwallJo{'2'^(^holeT^) 



(AlO) 

where Ohoie is the angular radius of the hole and Iwaii is the intensity of emission from the wall. The position of the 
null in the visibility function of a thin wall will then be 



'R'nuu{k\) 



77916 f D 



^source 



1 AU 



100 pc J \Rho 



(All) 



A. 3. Application to TW Hya 

The TW Hya disk has an outer radiu s of 196 AU (CO emission, IQi et all ((200l ') and an inner hole with radius 
4 AU (SED models, iCalvet et al.l ((200l : imaging, this paper). These size scales correspond to (u,v) distances of 
54 kA < 7?. < 2600 kA, a range well m atched to the coverage of the Very Large Array 7 millimeter observations, 
and we may apply the method of i jA.ll to generate the visibility function based on the simple power-law disk model. 
Approximating the Calvct ct al. (2002) model with power laws in surface density and temperature yields the profile 
T{R) = Two (R/ 100 AUy^ and T,{R) ^ J^^a {R/lOO AUy where Tioo= 28 K, g==0.44, I]ioo=3.7g cm^^ and 
p—0.90 For the dust opacity, we adopt a power- law distribution with = KQ{v/vf))^, where kq = 1.8cto^ 
vo — 1.0 X lO^Hz, and (3 — 0.8 ()D'Alessio et al.l 2 001 ). F igure [4| shows the visibility function calculated according 
to equation IA2I with the limits as in equations IA5I and IA7[ expanding to third order in (TZda)^- The resulting curve 
agrees well with the visibility function derived for the power law disk model obtained with the full radiative transfer 
calculation, shown in Figure [3 (light solid line) . Figure 2] also shows the visibility function of a thin ring interior to 
the disk (dotted line), as in ^A.2\ with radius 4 AU and flux 1.7 mJy. 
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